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• This talk gives a brief discussion of the motivations and principles behind computer 

04 ' automated expert semi-analytical resummation (caesar) for QCD final states. 

Global properties of energy- momentum flow in final-states, such as event shapes 
and jet-resolution thresholds, offer a good compromise between simplicity and sen- 
C*~) • sitivity to the dynamics of QCD. Thanks to the former, it has been possible to 

make a wide range of theoretical predictions for them, including both fixed-order 
and resummed perturbative calculations as well as non-perturbative model calcu- 
' lations. This set of theoretical tools has been essential in helping us to learn about 

. QCD from the extensive experimental data, having led in particular to numerous 

C*~) ' of measurements of the running coupling, a s pp , tests of the colour factors of QCD 

[2] and highly instructive studies concerning hadronisation (for reviews, see |3I4| ) . 
Most of the investigations so far have been carried out for observables that are 
. sensitive to the deviation from two-jet (1 + 1-jet) structure in e + e~ (DIS). This fairly 

straightforward environment has been of considerable value, especially in developing 
and refining our ideas about hadronisation. However the most stringent tests of 
the understanding thus gained, as well as possible novel extensions, are likely to be 
C^, found in the context of multi-jet events, including those at hadron-hadron colliders. 

A difficulty in such studies is that the next-to-leading logarithmic (NLL) re- 
summed calculations, needed in the limit of only soft and collinear emissions, have 
up to now usually been carried out by hand, a tedious and error-prone procedure, 
which has to be repeated separately for each new observable in each process. The 
' complexity of such calculations can increase substantially as one goes to multi-jet 

processes 

This is to be contrasted with fixed-order calculations, usually embodied in the 
form of a 'fixed-order Monte Carlo' program (FOMC, e.g. 6 ), which given a sub- 
routine for the observable, provides, numerically, the leading and next-to-leading 
order (LO, NLO) coefficients of the perturbative series. Though the internal com- 
plexity of FOMCs increases for processes with extra legs, the principles of their use 
remain identical regardless of the process. This makes it possible, even for those 
not expert in higher order QCD calculations, to obtain fixed-order predictions for 
arbitrary observables, in a wide range of processes. 

It is natural therefore to ask whether analogous programs could be constructed 
for resummed calculations. To help explain a solution to this question, proposed 
in 0, it is useful to define the problem a little more precisely. We will con- 
sider a general event shape (or jet-resolution threshold), defined by some function 
V(qi, . . . , gjv) of the momenta qi, . . . ,qrj of the N particles in an event. If the 
observable is such that it vanishes smoothly in the limit of n narrow jets, then we 
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call it an (n + l)-jet observable, and resummation will be needed in the n-jet limit. 
The resummed calculations should provide results for the probability f(v) that the 
event shape has a value smaller than v. For small v, this probability can often 
(exceptions are discussed below) be written in the 'exponentiated' form 



where Lg 1 (a s L) contains LL terms, a"L n+1 , g 2 {a s L) the NLL ones and so forth. 
Usually NLL accuracy represents the state of the art (though a related kind of 
resummation has recently been performed to NNLL accuracy {o, s g 3 ) for the e + e~ 
energy-energy correlation 8 ). 

In analogy with FOMCs, which calculate the LO and NLO contributions to f(v), 
the aim of an automated resummation program should be the calculation (sepa- 
rately) of the resummed LL and NLL contributions, i.e. the gi(a s L) and g2(a s L) 
functions, free of any contamination from potentially spurious higher orders, and 
in a form such that they can be expanded, allowing matching to fixed-order calcu- 
lations. 11 

The construction of an automated resummer, as least in our approach [7j, pro- 
ceeds however in quite a different fashion from an FOMC (usually just a Monte 
Carlo integrator with appropriate weights). This is because, in order to separate 
out the structure of Q it is useful to have some analytical understanding of the 
observable — for example regarding its behaviour in the presence of a soft and 
collinear emission, nearly always of the form 



where {p} are the hard momenta in the event and k is a soft emission that is 
collinear to the hard momentum (leg) with index £, with a transverse momentum, 
rapidity and azimuthal angle k t , r\ and <fi relative to leg £. For a large class of 
observables, given the coefficients ai and be, the function gi{a s L) has a known 
simple analytical form. So an automated resummer can determine the LL terms by 
verifying the parametrisation (J2J) and determining the coefficients. 

This sounds remarkably simple, and indeed there is a 'catch': one has also to 
establish that a given observable belongs to the class for which the LL terms truly 
are just determined by the values of the at, bg. Understanding precisely what that 
class is has been one of the significant developments of 0, with the introduction of 
a concept that we have called recursive infrared and collinear (rlRC) safety. Normal 
infrared and collinear safety (IRC) states that when an emission is very soft and/or 
collinear it should have no effect on the observable. The recursive variant essentially 
extends this condition to two-scale problems (as is relevant when resumming) and 
states that given an ensemble of arbitrarily soft and collinear emissions, the addition 

°This can be contrasted with the output from exclusive Monte Carlo event generators like Herwig 
151 or Pythia 1101 . which embody resummation for almost any observable via parton branching 
algorithms, and which provide an overall result for f(v), generally correct to LL accuracy, some- 
times to NLL, but for which it is difficult to guarantee the accuracy, extract separately the LL 
and NLL terms, perform matching to fixed order, and unambiguously separate the pcrturbativc 
and non-perturbative effects. 



f(v) ~ exp [Lgx{a s L) + g 2 (a s L) + a s g 3 (a s L) H 



L = ln~, (1) 
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of further emissions of similar softness or collinearity should not change the value 
of the observable by more than a factor of order one. Furthermore the addition of 
relatively much softer or more collinear emissions (whether with respect to the hard 
leg or one of the other emissions) should not change the value of the observable by 
more than some power of the relative extra softness or collinearity. While at first 
sight this seems similar to normal IRC safety, a number of observables have been 
found to be IRC safe but rlRC unsafe. Recursive IRC safety is a sufficient condition 
for the form Q to hold. While in principle one could also imagine resumming rlRC 
unsafe observables, we are not aware of any cases where this has been done. 

At NLL accuracy, it is not sufficient to know just the observable's behaviour 
when there is a single soft and collinear emission — additionally one needs to know 
how the observable behaves when there are multiple soft and collinear emissions 
of roughly similar hardness (its 'multiple-emission dependence'). Analytically this 
usually requires an in-depth analysis of the observable. However, as was shown in 
there exists a general solution to the problem, namely that at NLL all multiple- 
emission dependence can be taken into account via a contribution — \nT(R') to g2, 



where R'(a s L) = —dLLgi(a s L) (R' + 1 determines the average number of 'rele- 
vant' emissions), R' e is the part of R' associated with leg £, and Ki(v) is any soft 
momentum collinear to leg li that satisfies V({p},Ki(v)) — v. Since J-(R') can be 
calculated by Monte Carlo integration, one eliminates the need for any analytical 
study of the observable (the various limits can be taken numerically), thus providing 
the main NLL element of an automated resummation. 

One unsatisfactory issue that remains with the above approach is connected 
with the problem of globalness. It was shown in ^5] that for observables that are 
sensitive to emissions only in a restricted angular region (for example a single jet), 
additional classes of NLL terms arise, termed non-global logarithms (NGLs), so far 
calculated only in the large- N c limit (though see also J3j). Such NGLs depend 
non-trivially on the details of the boundaries of the restricted region in which the 
observable is sensitive to emissions, and these boundaries can be quite complex. 
Accordingly, for the time being, we have chosen to consider only the simpler case 
of (continuously) global observables, which have the characteristic that the de are 
all non-zero and the ae are all equal. These observables' resummations are free of 
NGLs. 

The resulting resummation program is called CAESAR, the Computer Automated 
Expert Semi- Analytical Resummer, and many results obtained with it can be found 
at http://qcd-caesar.org, for e + e~, DIS and hadron-hadron collisions. 

Particularly topical are hadron-collider dijet event-shapes, the subject also of 
recent experimental ^2] an d theoretical fixed-order ^5] work. These are of interest 
for a variety of reasons, ranging from the possibility of 'standard' measurements 
such as a s , to more novel studies examining for example the underlying event with 
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techniques such as those discussed in \IQ > or the soft (colour evolution) anomalous 
dimension matrices of ^7|. The event shape studied in |14I15| (a transverse thrust) 
has the drawback that its resummation involves NGLs, and so is beyond the current 
scope of CAESAR. To allow more straightforward compatibility with resummcd 
calculations, it is instead preferable to specifically design a set of global dijet event 
shapes (and jet rates). This has been the subject of work in O ne °f the main 
practical issues discussed there is the question of how to reconcile the requirement 
of globalness with the limited experimental detector coverage at large rapidities. 
Fortunately this problem can be circumvented in a variety of ways. This opens 
the road at hadron-colliders to full resummed studies matched with fixed order 
calculations [2j , with the prospect of significantly extending the very fruitful studies 
carried out in e + e~ and DIS over the past years. 
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